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Dark matter (DM) annihilation during hydrogen recombination {z ~ 1000) will alter the re¬ 
combination history of the Universe, and affect the observed CMB temperature and polarization 
fluctuations. Unlike other astrophysical probes of DM, this is free of the signihcant uncertain¬ 
ties in modelling galactic physics, and provides a method to detect and constrain the cosmological 
abundances of these particles. We parametrize the effect of DM annihilation as an injection of 
ionizing energy at a rate Cdm, and argue that this simple “on the spot” modihcation is a good 
approximation to the complicated interaction of the annihilation products with the photon-electron 
plasma. Generic models of DM do not change the redshift of recombination, but change the resid¬ 
ual ionization after recombination. This broadens the surface of last scattering, suppressing the 
temperature fluctuations and enhancing the polarization fluctuations. We use the temperature and 
polarization angular power spectra to measure these deviations from the standard recombination 
history, and therefore, indirectly probe DM annihilation. The modifications to the temperature 
power spectrum are nearly degenerate with the primordial scalar spectral index and amplitude; 
current CMB data are therefore unable to put any constraints on the annihilation power. This 
degeneracy is broken by polarization; Planck will have the sensitivity to measure annihilation power 
edm{z = 1000) > 10“^® eV/s/proton, while high sensitivity experiments (eg. NASA’s CMBPOL) 
could improve that constraint to edm{z = 1000) > 4 x 10“^® eV/s/proton, assuming a fractional 
detector sensitivity of AT/T ~ 1/rK and a beam of 3b These limits translate into a lower bound 
on the mass of the DM particle, Mdm > 10 — 100 GeV, assuming a single species with a cross 
section of {oav) ~ 2 x 10“^® cm®/s, and a fraction / ~ 0.1 — 1 of the rest mass energy used 
for ionization. The bounds for the WMAP 4y data are signihcantly lower, because of its lack of 
high S/N polarization measurements, but it can strongly constrain 0{ MeV) particles such as those 
proposed by Boehm et al (2004). 

PACS numbers: 


I. INTRODUCTION 

There is a broad consensus that the majority of the 
matter in the universe is non-luminous, non-baryonic 
“dark matter” (DM) [see Hi for recent reviews]. Some 
of the first compelling evidence for DM came from galaxy 
rotation curves, suggesting that a large fraction of the 
mass lay beyond the luminous extent of the galaxy. Simi¬ 
lar conclusions were reached for massive clusters of galax¬ 
ies using the gravitational distortion of the images of 
background galaxies. Measurements of the deuterium 
abundance, combined with big-bang nucleosynthesis in¬ 
dicate that density of baryons is less than the estimated 
mass density. This is supported by measurements of the 
temperature fluctuations of the cosmic microwave back¬ 
ground (CMB) radiation made by the WMAP satellitel^. 
Dark matter is an integral part of the cosmological “stan¬ 
dard model”, constituting 80% of the total matter in the 
Universe. 

Cosmology also constrains the gravitational properties 
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of the DM. The spatial clustering of galaxies and the 
angular power spectrum of CMB temperature fluctua¬ 
tions strongly favor a non-relativistic pressureless species 
for the DM, while N-body simulations show that this 
clustering is consistent with the DM being composed of 
weakly interacting particles whose dominant long-range 
interactions are gravitational. This suggests that the DM 
particle is a weakly interacting massive particle (WIMP), 
although more exotic possibilities are not ruled out. We, 
however, will focus solely on WIMPs for this paper. 

A theoretical description of the DM is still unknown, 
although it is believed that such a description involves 
physics beyond the Standard Model. There are strong 
theoretical reasons for believing that the Standard Model 
is modified at the electroweak symmetry scale ^ 1 TeV, 
intriguingly the same energy scale for WIMPs predicted 
by the present mass density. These modifications range 
from supersymmetric extensions of the Standard Model 
to large extra dimensions modifying gravity on these 
scales, and generically have a zoo of massive particles 
with properties that make them potential DM candi¬ 
dates. The relevant energies are just entering the reach 
of current direct detection experiments^ and accelera¬ 
tor energies, and will be strongly constrained by the next 
generation of these experiments. 

While there is no substitute for a direct detection of 
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WIMPs, astrophysical probes of DM play an important 
role, since they provide a complementary view of the 
DM parameter space, making different assumptions than 
particle physics experiments. Furthermore, probing the 
cosmological abundance of a candidate particle requires 
an astrophysical probe. There already are a few tan¬ 
talizing observations suggesting a DM particle between 
~ 1 — 100 GeV. The 7 -ray emission measured by the 
EGRET instrument on the Compton Gamma Ray Ob¬ 
servatory has a higher amplitude at ^1 GeV than tra¬ 
ditional models of the Galactic cosmic rw population 
and interstellar medium models allow S 0 ■ The 
annihilation line strength at 511 keV observed by Inte¬ 
gral/SPI suggests a surprising Galaxy-wide positron pro¬ 
duction rate of 10“^^ s“^ 0. Furthermore, the cosmic- 
ray positron ratio observed by HEAT shows an excess 
above 5 GeV consistent with simple models of Higgsino 
decay 0. The synchrotron haze in the inner Milky Way 
|l 0 | may be an example of synchrotron emission from 
pairs produced by ongoing DM annihilations as sug¬ 
gested by M The observed synchrotron signal could be 
produced by a hducial model of 100 GeV particles anni¬ 
hilating at (aAv) = 2 X 10“^® cm^ s“^ and distributed 
with an NFWjl^ mass profile [T^. This model is not 
expected to be correct in detail, but is a hducial model 
scalable to other possible scenarios. 

In this paper, we propose using the GMB tempera¬ 
ture and polarization fluctuations as a probe of DM an¬ 
nihilation. DM annihilation at 2 ; ~ 1000 injects energy 
into the photon-baryon plasma ionizing neutral hydrogen 
and modihes the recombination history. These additional 
electrons scatter GMB photons, making the last scatter¬ 
ing surface thicker and attenuating correlations between 
temperature perturbations. On the other hand, the cor¬ 
relations between polarization fluctuations are enhanced 
by the thicker scattering surface. This alters the temper¬ 
ature and polarization angular power spectra, providing a 
handle on the properties of the DM. The GMB has an im¬ 
portant advantage over the other probes discussed above, 
as computing the power spectrum is a linear calculation 
based on well understood physics. This enables detect¬ 
ing small deviations from the expected signal (that herald 
new physics) at a high level of signihcance. This should 
be contrasted with the other probes discussed above that 
require additional information such as the DM distribu¬ 
tion and dumpiness, ISM density, magnetic field strength 
and degree of tangling, Galactic photon energy density 
etc., all of which are complex processes with significant 
uncertainties. 

Modifications to the recombination history, ranging 
from delayed recombination to low redshift ionization 
from DM decays, and their effect of the CMB have been 
considered by [eg. El III 13111111113 Most recently, 

these ideas have been used to attempt to explain the 
high optical depth observed by the WMAP satellite by 
using DM decays to reionize the universe. The goals of 
this paper significantly differ from such studies. We aim 
to understand how the GMB (via the recombination his¬ 


tory) constrains the parameters of an annihilating DM 
model. Since we lack a preferred theoretical framework 
for the DM, we construct a generic parametrization of 
DM annihilation to effectively constrain the DM model 
space. We parametrize DM annihilation as an energy 
injection into the IGM^l and relate this energy to stan¬ 
dard parameters like the DM mass and annihilation cross 
section. Given a specific model of the DM, constraints 
on this injection energy can then be used to constrain the 
parameters of that model. 

We start (Sec. InJ by discussing the effect of DM an¬ 
nihilation on the ionization history of the Universe, in¬ 
troducing the “on the spot” approximation as a useful 
parametrization of the more complex processes responsi¬ 
ble. Sec. lIIIl then discusses the effect of this change to the 
ionization history on the GMB, while Sec. IIVI attempts to 
quantify the detectability of these effects. Sec. m consid¬ 
ers whether other astrophysical probes have the poten¬ 
tial to probe DM annihilation during the recombination 
epoch, and Sec. I VII summarizes the principal conclusions 
of this paper, and discusses their implications. 

A note on cosmology and notation: if not explic¬ 
itly specified, we will assume a concordance fiducial 
model with + Da = 1, Dm = 0.3, Df, = 0.05, 
h — 0.7, Yffe = 0.24 and ris = 1. We only con¬ 
sider the scalar contributions to the GMB fluctua¬ 
tions, and therefore describe the temperature and po¬ 
larization fluctuations by the temperature-temperature 
(TT), temperature-polarization (TE), and polarization- 
polarization (EE) power spectra. Finally, proper and co¬ 
moving times are t and ry respectively, while a 0 subscript 
denotes the present epoch. 


II. THE IONIZATION HISTORY 

The effect of DM annihilation on the recombination 
history can be conveniently, albeit artificially, separated 
into two stages - the injection of the energy from the an¬ 
nihilation into the IGM, and the effect of this energy on 
recombination. We argue that the former process is well 
approximated by a rate of energy injection per hydro¬ 
gen atom (edm), instantaneously used to ionize and heat 
the IGM. We then consider how this energy changes the 
ionization fraction as a function of time. 


A. The “On the Spot” Approximation 

How do DM annihilations cause ionizations? The pri¬ 
mary products from an annihilation depend on the par¬ 
ticular DM model, but generically are quarks, gauge 
bosons, leptons, and Higgs particles. These primaries 
tend to be unstable, and rapidly decay via hadronic- 
leptonic jets into showers of pairs, protons, photons 
and neutrinos. Given our ignorance about DM, we as¬ 
sume each annihilation partitions the majority of its en¬ 
ergy between pairs, photons and neutrinos, whose en- 
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FIG. 1: A comparison of the photon cooling time to the 
Hubble time at 2 = 1000, for different photon energies. The 
dominant processes (in order of increasing energy) are ion¬ 
ization, Compton scattering, pair production, and photon- 
photon scattering. All the curves (except the dotted curve) 
assume a neutral IGM, with a density of 2 x 10”'" cm“® atoms 
today. The dotted curve shows the pair production rate for 
a completely ionized IGM. Regions where tH/tcooi < 1 are 
transparent; photons injected at these energies lose their en¬ 
ergy by redshifting. Note that photon-photon scattering does 
not transfer the energy to electrons, but simply redistributes 
it to lower energies. This hgure ignores pair production off 
CMB photons since this process is subdominant for the en¬ 
ergy range considered here; it however dominates at higher 
energies. 


ergy spectra can be calculated given the mass and cou¬ 
plings of the DM particle. Note that we are not assum¬ 
ing that these are directly produced by the annihilation, 
but simply that they are final products of the resulting 
particle cascades. The problem now simplifies to under¬ 
standing the mechanisms by which pairs, photons and 
neutrinos inject energy into the IGM. Of these, neutri¬ 
nos are the easiest to understand; they never interact and 
their energy is lost. 

The interaction of photons with the IGM was consid¬ 
ered in detail by pn| who find that the dominant pro¬ 
cesses (ordered by increasing photon energy) are pho¬ 
toionization, Gompton scattering, pair production off nu¬ 
clei and atoms, photon-photon scattering, and pair pro¬ 
duction off CMB photons (Fig. ^|. To estimate the ef¬ 
ficiency of these mechanisms, we compare the cooling 
time for each process, tcooi = \/{d\nE/dt), to the Hub¬ 
ble time, tn = 1/H{z). Except for Gompton scattering, 
we approximate the cooling time by the mean free time 
as most of the energy is lost in the first interaction. If 
tn > tcooi (Fig. d, energy deposition is very efficient, 
either by directly ionizing the IGM, or by producing en¬ 
ergetic electrons. Conversely, if tn ^ tcooi, the universe 


is optically thin and most of the energy is lost through 
the redshifting of photons and not to ionizations. The 
photon-photon scattering process is an exception to 
the above - each scattering event, on average, equally di¬ 
vides the energy between the two photons. The photon 
energy spectrum therefore gets shifted to lower energies 
until either pair production starts to dominate, or the 
universe becomes transparent. 

What happens to photons injected into the trans¬ 
parency window between ~ 10 ® — GeV (Fig. d? 
The ratio tu/tcooi oc (1 -I- (cx (1 -I- 2 )®/^ for two 
photon scattering), while the photon energy redshifts as 
(1 -I- 2 ). These photons therefore remain in the optically 
thin regime, and contribute to the diffuse photon back¬ 
ground today (Sec. 

The second component of energy injection comes from 
electrons ^ 3 , both from the annihilation products, as 
well as from Compton scattering and pair production 
considered above. Their energy loss has been considered 
by a number of authors [ij, ; we restrict ourselves 

to a brief discussion of the relevant processes and time 
scales. At high electron energies (7 1 ), the dominant 

energy loss is by inverse Compton scattering CMB pho¬ 
tons. The cooling time is |^ . 

1 ^ -d\n-i ^ 

tcooi J dt ^TTlcC^ 


where Tcmb = 2.725(1-|-2)K is the mean CMB tempera¬ 
ture at the relevant redshift, an is the radiation constant, 
and tJT = 6.65 x cm^ is the Thomson cross section. 
Comparing this to the Hubble time, one finds 


tn 


tcooi 


10 ® 


1 + 2 ^/^ 1 
1000 ) ’ 


( 2 ) 


implying that inverse Compton cooling efficiently pro¬ 
duces photons with energies 


Ej 



Ec 


1 GeV 


2 

MeV . 


( 3 ) 


Fig.dshows that electrons with energies < 100 MeV will 
produce photons that efficiently ionize hydrogen; above 
that energy, the scattered photons either produce an elec¬ 
tromagnetic cascade by Compton scattering or pair pro¬ 
duction, or are scattered into the optically thin part of 
the spectrum from 10® to 10^^° eV and escape. 

At lower energies, the principal mechanisms for en¬ 
ergy loss become collisional heating, excitations and 
ionizationsp^. At high kinetic energies, Eg 100 eV, 
the cross section for collisional ionization is 


2.23 X 10-1® ln(F;/13.6) , 

cTeH = - ^ - -cm^ , (4) 

h/ 

{E measured in eV) implying tn/tcooi ^ 1 at 2 ~ 1000. 
The results for collisions and excitations are very simi¬ 
lar, although collisional losses become important at lower 
energies, 100 eV < if < 1 keV. 
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FIG. 2: The injection of energy from dark matter annihilation 
into the IGM, via the creation of electromagnetic cascades. 
Energy transfer to the IGM takes place principally through 
the ionization and collisional processes. 


Our model of DM annihilation, summarized in Fig. |21 

1. DM particles annihilate into jets, whose end prod¬ 
ucts are dominated by electrons, photons and neu¬ 
trinos. 

2. The IGM is transparent to neutrinos; this energy 
is lost. 


3. The photons and electrons trigger electromagnetic 
cascades, shifting their spectra to lower energies, 
until their energy is either deposited in the IGM, 
or is redshifted away when the photons enter an 
optically thin regime. 

4. The time scales for the cascades and energy depo¬ 
sition are much smaller than the expansion time. 

We parametrize the effect of DM annihilation by the rate 
of energy injection per hydrogen nucleus per time, edm', 
furthermore, we assume this energy is instantaneously 
deposited into the IGM. This “on the spot” approxima¬ 
tion has the virtue of being generic and independent of 
particular properties of the DM. 

What is the magnitude and redshift dependence of 
Cdm? Given a particle with mass Mdm and a thermally 
averaged cross-section {(tav), we obtain. 


Crfm — f fbfrfn 




'dm^Q 




(l + z)3 


( 5 ) 


where noMfi and nnfi are the present densities of the 
dark matter and hydrogen particles respectively, and / 
is the fraction of the rest mass energy injected into the 
IGM. Assuming our fiducial cosmology with a single DM 
species, we find. 


f-dm ~ -k zf eV s~^ 

lOOGeVW (aAv) V 

Mdm / \2 X cm3 


( 6 ) 


The injected energy is inversely proportional to the par¬ 
ticle mass; more massive particles inject less energy into 
the IGM. We parametrize our ignorance of the annihila¬ 
tions and their effect on the IGM by a simple efficiency 
factor, /. Given a specific model, one can compute / and 
convert constraints on edm into constraints on Mdm and 
other model parameters. 


B. Recombination with DM annihilation 


Given Cdm, we compute its effect on the recombina¬ 
tion history. This energy injection heats the IGM, and 
ionizes and excites the hydrogen and helium atoms. 
compute the exact fractions converted to heat, ionization 
and excitation as a function of the ionization fraction, 
and find that for a neutral IGM, the energy is roughly 
equipartitioned between the three processes, while for a 
fully ionized plasma, all the energy is converted into heat. 
This suggests a simpler but adequate approximation M 
that (1 — x)/3 of the energy goes into ionization and 
(l-|-2a;)/3 into heating the IGM, where x is the ionization 
fraction and we assume that excitations neither change 
the matter temperature nor the ionization fraction. 

We compute the recombination history using the pub¬ 
lic code RECFAST modifying the evolution equations 
as follows. 





^dm,0 1 \ 

13.6 3(1 + fHe) 


^dm,o 1 . 

24.6 3{l + fHe) 


( 7 ) 


where edm,o is the energy injection rate at the present 
epoch (in eV/s), and 


^ (1 + ^) 
H{z){l + z) • 


Additionally, the ionization fraction of a species A is de¬ 
fined as. 


a:[A] = 


n[A+] 

n[A+] -k n[A\ ’ 


( 9 ) 


and fHe is the ratio of the number density of helium to 
that of hydrogen. The evolution of the matter tempera¬ 
ture, Tm is similarly given by, 


JdTm\ 

[nr) = 

l + 2,T[fl| + /H,(l + 2i|He])„ . 

-sTs - wns - 


X 
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FIG. 3: The ionization fraction Xe (top), matter tempera¬ 
ture (center), and visibility function (bottom) as a function 
of Cdm- The heavy solid lines show the fiducial model with 
Edm = 0; from bottom to top, Cdm.o = 5,10, 100, 500 x 10“^® 
eV/s. The thin dashed line in the center plot shows the evo¬ 
lution of CMB temperature, T{z) = ro(l -I- z). Note that the 
injection of additional energy does not slow recombination, 
but increases the residual ionization; this leaves the peak of 
the visibility function unchanged but broadens the surface of 
last scattering. 


The resulting recombination and matter temperature 
histories for different values of edm,o are shown in Fig. 0 
the dominant effect is to change the residual ionization 
after recombination. This is easily explained by consider¬ 
ing the competition between the recombination rate and 
the expansion of the universe. At early times, the recom¬ 
bination rate is significantly greater than the expansion 
rate and therefore, additional ionizations due to DM an¬ 
nihilation are immediately erased. As the recombination 
rate slows, these additional ionizations “freeze out”, lead¬ 
ing to a greater residual ionization fraction. 

The evolution of the matter temperature is similar. 
At redshifts ^ 100, Compton scattering keeps the mat¬ 
ter and radiation in tight thermal contact, and the ex¬ 
cess energy from DM annihilation is lost in the extremely 
large heat capacity of the blackbody radiation. However, 
as the matter completely decouples from the radiation, 
annihilations start to increase the matter temperature, 
resulting in slower cooling relative to the fiducial model. 


III. THE CMB AS A PROBE 

Having computed the effect of DM annihilation on the 
recombination history, we attempt to understand its ef¬ 
fect on the CMB. In what follows, it is sometimes con¬ 
venient to parametrize the effect of DM annihilation by 
an ionization “floor” added to the standard recombina¬ 
tion history. This lacks the physical intuition of Cdm, but 
is a convenient analytic approximation. We define the 
optical depth to Thomson scattering, 

rvo 

T{r]) = / dijaTrieCa , ( 11 ) 

Jr) 

where Ue is the free electron density. Assuming a mat¬ 
ter dominated cosmology and constant ionization frac¬ 
tion Xe, this gives us 

t { z ) ^ 4 X ^ (^ 3 ) 

if 1. 

A. Peak Positions 

We begin by estimating the change in the position of 
the acoustic peaks in the temperature power spectrum 
due to an ionization floor. The probability that a photon 
last scattered between redshifts 0 and z dz is given by 
the visibility function, 

g{z) = T'(z)e-"(^) , (13) 

shown in Fig. |21 for different recombination histories. 
The fraction of photons that scatter at a redshift < z, 
G{z), is simply the integral of visibility function, G{z) = 
1 — exp(—r(z)). Since g{z) is sharply peaked, we can 
meaningfully define a redshift of last scattering, zls , that 
determines the angular positions of the acoustic peaks. 
A convenient definition is G{zls) = 0-5 or t{zls) ~ 0-7 
implying zls ~ 1050 for standard recombination. For 
the ionization floor to significantly shift the peaks, the 
additional optical depth, At, would have to be ~ 1. Us¬ 
ing Eq. 113 this requires 10^Xejioor ~ {V^m/^ bh), or 
Xe, floor ~ 0.01 for our fiducial cosmology. As we shall 
see below, such an ionization fraction would have already 
noticeably affected the CMB temperature and polariza¬ 
tion and therefore is strongly disfavored. More plausible 
values of the ionization floor do not noticeably shift the 
positions of the acoustic peaks in the temperature power 
spectrum. 

B. Power Spectra 

The effect of the altered recombination history on 
the CMB power spectra is discussed analytically below. 
However, the numerical results presented in the paper use 
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FIG. 4: The TT, TE, and EE angular power spectra for our 
fiducial cosmological model, with no DM annihilation (solid 
and dotted lines), and with edm,o = 10“^^ eV/s. Also shown 
are the polarization noise spectra, for the WMAP V band, the 
Planck 143 Ghz channel, and a hypothetical high resolution 
polarization experiment (see Tabled] for details). 


the publicly available Boltzmann code CAMB , with the 
modified version of RECFAST described in the previous 
section, to obtain accurate power spectra. An example 
is shown in Fig. 0| 

The temperature angular power spectrum is the pho¬ 
ton distribution function convolved with the visibility 
function (the last scattering surface), and projected on 
the sky. The photon distribution function is unchanged 
by DM annihilation, but the visibility function extends to 
lower redshifts, broadening the surface of last scattering. 
This suppresses perturbations on scales smaller than the 
width of the surface, resulting in a relative attenuation 
of the power spectrum. This is scale dependent, with 
the largest scales attenuated the least and small scales 
the most. These effects are clearly seen in the accurate 
numerical solutions in Fig. 21 

Given the imminent high S/N temperature measure¬ 
ments due from the WMAP and Planck satellites, an 
immediate question is whether DM annihilation is de¬ 
tectable just using the temperature power spectrum. Un¬ 
fortunately, the effects of e^m described above are almost 
perfectly degenerate with the slope and amplitude of the 
primordial power spectrum. To see this quantitatively, 
we start with the line of sight solution to the tempera¬ 
ture perturbation in direction n [^ . 


A(fi) 



T' 


0 


n • Vf, 


-l- 20 


e dr] , (14) 


where Ih is the gravitational potential, 0 is the photon 
density perturbation, V{, is the baryon velocity, and we 
ignore vector and tensor contributions. If we ignore the 


ISW [23 contribution (20) to the anisotropy spectrum, 
we obtain a useful semi-analytic approximation to the 
anisotropy spectrum by separating into slowly varying 
(potentials, T{k) below) and rapidly var ying (recombi¬ 
nation, Silk damping, D{k) below) terms |2g, 

nOO 

Ci=AttA d{lnk)k^‘D^{k)T^{k) , (15) 

Jo 

implicitly assuming that T‘^{k) is evaluated at the red- 
shift of last scattering and has no time dependence. The 
damping function is given by [^ . 

D{k)=J <iz5(z)exp , (16) 


where g{z) is the visibility function introduced earlier, 
and kn is the Silk damping scale given by 

j_ _ A°° _c_ 1 r 

kjj ^ H'^{z) 6{1 + R)t'{z) (1-I-i?) 15 

(17) 

where R = 3pi,/4p^ is the baryon-photon ratio. Since 
the ionization history only appears in Ea. ll5l through the 
optical depth in D{k), we estimate the effect of adding an 
ionization floor by computing D{k)/DQ{k), where D^ik) 
assumes the standard ionization history. 

As the relevant regime is when the ionization fraction 
is rapidly changing, we numerically integrate Eq. d and 
compute D{k)/Do{k) for different edm,o- The results for 
our fiducial cosmology are shown in Fig. 21 The scales rel¬ 
evant for Z > 50 in the CMB correspond approximately to 
k > O.OOl/iMpc”^; Fig.Odemonstrates that D{k)/Da{k) 
is remarkably well described by a power law, over 
these scales. This signals a near exact degeneracy in the 
CMB; examining Eq. suggests that the effect of the 
ionization floor can be almost exactly compensated by 
adjusting ris ^ ris + 2a, and changing the amplitude, A. 
The residual differences can be corrected by adjusting (at 
sub-percent levels) the remaining cosmological parame¬ 
ters. 

We emphasize that this degeneracy appears to be 
purely accidental. As fc ^ 0, Silk damping becomes in¬ 
creasingly unimportant and D{k)/Do(k) 1. In addi¬ 
tion, we have ignored the ISW contribution, which has 
a different visibility function, and therefore will not be 
compensated by changing the scalar spectral index. On 
small scales {k —> 00 ) that are considerably damped be¬ 
fore recombination, the correction to the visibility func¬ 
tion due to the ionization floor is negligible and again, 
one would expect D{k)/Do{k) ~ constant. These two 
limits are however hard to constrain, large scales be¬ 
cause of cosmic variance, and small scales because of sec¬ 
ondary anisotropies. On intermediate scales where high 
S/N measurements of CMB can be made, the effect of 
an ionization floor is degenerate with changing the scalar 
spectral index. 

Estimating the effect of Cdm on the polarization of the 
CMB is more involved. Polarization principally results 
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the monopole, resulting in the well known phase struc¬ 
ture of the CMB temperature and polarization spectrum 
peaks. As the last scattering surface grows thicker, the 
fractional contribution from monopole perturbations to 
the quadrupole increases, shifting the positions of the TE 
and EE peaks. Finally, on smaller scales, the TE and 
EE power spectra are attenuated by increased scatter¬ 
ing, analogous to the TT power spectrum. These trends 
are seen in the TE and EE power spectra in Fig.^ 



FIG. 5: (Top) The heavy lines show the ratio of damping 
functions D{k)/Do{k) for edm,o of 100, (solid), 500 (dotted), 
and 1000 x 10“^®eV/s (dashed) and onr fiducial cosmology. 
The light lines show this same ratio divided by {k/kfid)~°‘ 
where kfid = 0.05 Mpc“^. The horizontal dotted lines are 
a visual guide. (Bottom) The solid line shows the ratio of a 
model with no DM annihilation, but with Us altered using the 
analytic calculation above, to our fiducial model with edmfi = 
500 X 10“^®eV/s. The dashed line has the same ratio, except 
that all the cosmological parameters are adjusted to best fit 
the model with DM annihilation. 


from the Thomson scattering of the local quadrupole in 
the temperature distribution. However, the quadrupole 
vanishes during the tightly coupled regime before recom¬ 
bination; the only source of a quadrupole is the free 
streaming of the monopole and dipole perturbations dur¬ 
ing recombination. Ignoring the effects of reionization, 
the amplitude of the quadrupole contributing to polar¬ 
ization can be schematically written as |^ . 


IV. ESTIMATING DETECTABILITY 
A. Formalism 


Given a measurement of the CMB sky with detector 
noise and cosmic variance, how distinguishable are any 
two models? And does there exist a combination of stan¬ 
dard CMB parameters that can mimic DM annihilation? 
To answer the first question, we assume a realization 
of the full sky both in temperature and polarization (E 
modes). The likelihood of observing these maps is. 


/:(d|theory) oc-^=L=exp , ( 20 ) 

where we have assembled the maps into the vector d, and 
the covariance matrix, C, is a function of the theoretical 
model. We transform to the spherical harmonic basis. 


£(d I theory) oc 


n 


1 

Vdet Cl 


exp 



cr 



( 21 ) 

where = {aj^, af^) is the spherical transform of the 
temperature and E-mode maps respectively. The 2x2 
covariance matrix for each Z, m mode is given by 


Cl = 


cr 

cr 


cr 

qEE 


( 22 ) 


02 (fc) 


Qo{k)j2{x) + 301 (fc) 


ji{x) 


‘^h{x) 

X 


(18) 


where 0 ; with I = 0 , 1 , 2 represents the monopole, dipole 
and quadrupole components of photon distribution, and 
X = kAr], where Arj is the thickness of the last scat¬ 
tering surface. Focusing on scales much larger than the 
thickness of the last scattering surface, x <C 1 , we obtain. 


02(fc) 


eo{k)[kAr + 6Qi[kAr]] 3 

1 t; T ' ’ 


(19) 


where we used the expansion ji{x) = x’’/{2l + 1 )!! -|- 
0{x^^'^). Increasing the width of the last scattering 
surface therefore increases the amplitude of the polar¬ 
ization fluctuations. Furthermore, Eq. m implies that 
the quadrupole is dominated by free-streaming from the 
dipole perturbations. These are 7 r /2 out of phase with 


where C'®'®, and are the temperature and E- 
mode angular auto-power spectra and cross-power spec¬ 
trum respectively. Taking logarithms and summing over 
azimuthal modes, we get 


iog[/:(d|c)] = --^( 2 z + i) 


log(det C)-l- 


Ql 1 (j-bJ-bJ _|_ (jbdbbfjb b _ 

det C 


(23) 


where det C = — (C^^)^, the hat denotes ob¬ 

served quantities, and the I dependence is suppressed. 

We can now compare the likelihoods for two different 
theoretical models, C and C', 


r = log 


V/:(d|C0y' 


(24) 



























obtaining 


- -9 


log 


det C 


aC 


TT 


det C' 


where 


^=1 
7 =-2 


f^EE 

^EE' 

det C 

det a 

CTT 

CTT' ' 

det C 

det C' 

' 

CTE' 

det C 

det C' 


(25) 

(26) 

(27) 

(28) 


If we assume that the data is a realization of C, we have a 
simple criterion for distinguishability (and therefore, de¬ 
tectability) - two models are distinguishable if the proba¬ 
bility that r < 0, P(r < 0), is less than a chosen threshold 
(the confidence level). To compute P{r < 0), we compute 


1 


- -9 + 


log 


det C 
det C' 


-f -f 7C 


EE 


-,TE 


(29) 


and 


Var(r) = + 1) 




[(C^ ^ ] + 2af3{C^^y + 


+ + 2a'yC^^C 


TEr^TT 


2{r) , (30) 


where we use standard contraction formulae to compute 
the four point functions, and the last approximation is 
good when C is very close to C'. The central limit the¬ 
orem ensures that the distribution of r is well approxi¬ 
mated by a Gaussian; we therefore obtain. 


P(r < 0) = - 


1 — erf 


•\/ 2 Var(r) 


(31) 


The above expressions can be generalized to take into 
account of detector noise; one substitutes ^ + 

N^, where x = TT,EE assuming the temperature and 
polarization measurements are uncorrelated. The noise 
is determined by the angular size of the telescope beam, 
and the temperature sensitivity of the detectors pm, 

N{l) = {wp)-^exp[l{l + l)e^] , (32) 

where 9 is related to the FWHM of the beam by 
FWHM = 9V8h[2, (wp)"^/^ = AT x FWHM, and all 
angles are in radians. 

We now state our algorithm for the detectability of DM 
annihilation: 


Experiment Beam 10'’Ar/r 10'’Ar/r 

FWHM (arcmin) (I) (Q,U) 


WMAP (V band) 

21 

11.0 

15.6 

Planck (143 Ghz) 

7.1 

2.2 

4.2 

HiRes POL 

3.0 

1.0 

1.0 

Cosmic Variance 

0.0 

0.0 

0.0 


TABLE I: Detector sensitivities and beams for different CMB 
temperature and polarization experiments. HiRes POL refers 
to a hypothetical all sky CMB polarization experiment. 


Experiment 

68% 

edm,o(10"" eV/s) 
90% 

99% 

WMAP (V band) 1.5e2 

5.6e2 

l.le3 

Planck (143 Ghz) l.lel 

3.2el 

6.2el 

HiRes POL 4.1 

1.2el 

2.2el 

Cosmic Variance 4.0 

l.lel 

2.lei 


TABLE IL Detectable values of tdm at 68%, 90%, and 99% 
confidence levels, for different experimental parameters (Ta¬ 
ble EJ. 


1. Consider two cosmological models, one with Cdm = 
0 and the other with Cdm > 0. For simplicity, we 
assume a minimal cosmological model with 6 pa¬ 
rameters (Dm, Hf,, h, Us, A, r). 

2. We adjust these parameters for the model with no 
DM annihilation to minimize (r). 

3. At the minimum, we compute Var(r) and P(r < 0), 
thereby obtaining a measure of the detectability of 

^dm ■ 

This algorithm is complementary to the Fisher informa¬ 
tion methods currently popular in cosmology. The Fisher 
information probes the likelihood function in the neigh¬ 
borhood of a fiducial point, estimating the minimal the¬ 
oretically achievable errors. On the other hand, we ex¬ 
plicitly track the degeneracy locus well beyond the im¬ 
mediate neighborhood of a fiducial model. 


B. Results 

We consider four different experiments - a cosmic 
variance limited experiment, a hypothetical high reso¬ 
lution polarization experiment (hereafter HiRes POL), 
the Planck 143 Ghz channel and the WMAP V band 
Ttie parameters chosen for HiRes POL are sim¬ 
ilar to those projected for (eg. NASA’s CMBpol |4^1 
missions designed to detect the imprint of gravitational 
radiation from inflation. The assumed noise and beam 
characteristics of these experiments are in Table 13 and 
the noise power spectra are plotted in Fig. ^ Note that 
our aim here is to consider representative parameters, 
and not an optimization of experimental specifications. 

The results are summarized in Fig. El and Table HD 
A principal feature is the significant improvement of 
Planck over WMAP , due to the introduction of high S/N 
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FIG. 6: The detectability of a given injection energy, edm,o, 
for the different experimental specifications of TableEl assum¬ 
ing a maximum multipole, Imax = 2500. The lines, from left 
to right, are for cosmic variance, HiRes POL, the Planck 143 
Ghz channel, and the WMAP V band. The dashed line shows 
the cosmic variance detectability for Imax = 1500, while the 
dot-dashed lines show it for 20% and 50% sky coverage. The 
dotted line uses the parameters for HiRes POL, but with no 
polarization information. Also shown are the 68%, 90%, and 
99% levels. 


polarization data, evident in Fig. 21 The importance of 
polarization is further emphasized by comparing HiRes 
POL, to an experiment with the same temperature but no 
polarization sensitivity. As discussed earlier, the temper¬ 
ature power spectrum suffers from a degeneracy between 
Us and tdm, and therefore cannot constrain Cdm by itself. 
Introducing polarization breaks this degeneracy, allowing 
Cdm to be measured with significantly greater sensitivity. 
We also observe that HiRes POL (almost) achieves the 
cosmic variance sensitivity limits. 

The advantage of phrasing the limits in terms of edm,o 
is that they are independent of a particular DM model, 
or indeed, of any mechanism for the injection of addi¬ 
tional energy during recombination. However, for a sin¬ 
gle species of DM, Ep.^lrelates Mjjm, ^dm, {o'av) and / 
to tdmfi- Furthermore, / and (ctav) are degenerate with 
each other, allowing us to translate our limits on e^m into 
constraints in the J-Mum plane. Fig.Qdoes this for our 
fiducial cosmological model, using Eq. Eland Tableim We 
observe that, assuming / ~ 0.1 — 1, Planck will be able 
to detect (at 90% confidence) DM annihilation from par¬ 
ticles with masses less than ~ 3 — 30 GeV, while HiRes 
POL increases that lower bound to ^ 10 — 100 GeV. 
These limits assume (aAv) = 2 x 10“^® cm^/s, appro¬ 
priate for a thermal relic; a higher {uav) (due to eg. 
co-annihilations) would proportionally increase the lim¬ 
its. We note that these limits probe relevant parts of 


FIG. 7: The 90% exclusion region in the Mdm — f plane, 
for (from top to bottom) WMAP , Planck and HiRes POL; 
regions above the lines can be excluded by these experiments. 
This assumes our Hducial cosmology, and Eq. Elto relate / and 
Mdm to Cdm', (cav) is set to the value required for a thermal 
relic density, but deviations from this can be absorbed into 
/. The dotted lines also show the 68% and 99% exclusion 
regions for our hypothetical polarization experiment. 


parameter space, and are complementary to accelerator 
and direct detections, since they make very different as¬ 
sumptions. 

It is timely to ask what parts of DM model space can 
be constrained by WMAP . Fig.[7|shows that WMAP will 
have no sensitivity to DM models with masses > 1 GeV, 
and therefore, to the traditional DM candidates. How¬ 
ever, 113,0,113 have proposed a light (~ 10 — 100 MeV) 
DM particle to explain the 511 keV flux observed by the 
INTEGRAL satellite. Since the particle mass is below 
the photon transparency window, the annihilation en¬ 
ergy is efficiently converted into ionizations, suggesting 
that such models will be strongly constrained by the 
WMAP polarization measurements. 

The discussion above focuses on the detection of an 
unknown particle, and we therefore have concentrated 
on the simplest case of a single species of DM. Reality 
could well be more complicated, and the DM could well 
consist of multiple species. In such a case, in addition 
to determining the nature of the DM, we would need 
to determine the relative contribution of the different 
species. Given a particle with a known mass and an¬ 
nihilation cross section, constraints on Cdm put an upper 
bound on the density of these particles, without making 
any assumptions about whether the particle is a thermal 
relic or not. Eig. El demonstrates this for different values 
of Mdm, {o’av) and /. This emphasizes the complemen¬ 
tary nature of astrophysical and particle physics probes 
of DM; at low masses, the GMB is competitive with di- 
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FIG. 8: Contours in the Q,dm-{oav) plane with edm,o = 
10-=^ eV/s, for Mom// = 1 TeV, 1 GeV, and 1 MeV (solid), 
and 10 TeV, 10 GeV, and 10 MeV (dashed). Regions above 
these contours are accessible to an experiment with the sen¬ 
sitivity to measure edm,o = 10“^“^ eV/s. The shaded region 
shows the region excluded by our fiducial model. 


rect detection and accelerator searches of DM. However, 
given a detected DM particle, the CMB probes, with min¬ 
imal assumptions, its cosmological density, that in turn, 
has the potential to constrain theoretical models for its 
formation. 


V. OTHER ASTROPHYSICAL CONSTRAINTS 


Does DM annihilation during the recombination epoch 
have other observational consequences? We consider 
three possibilities - distortions to the CMB spectrum, 
redshifted photons contributing to the diffuse photon 
background, and molecular hydrogen production. 

We start by considering the additional energy density 
injected per unit time at z = 1000 due to DM annihila¬ 
tion. 


~ 2 X 10“^^ eV/cm^/s 


1 + zV 
1000 ) ^ 


0.12 


(ctav) 


100 GeV 


2 X 10 cm^/s Mdm 


Over a Hubble time, the energy injected is 


E - 2 eV/cm' 


^DAlh^ 


/ 0.147 


f—) f 

V 1000 J \ 


Hm/iV 


(ctav) 


.(33) 


100 GeV\ 


0.12 ) 1^2 X 10-26 cm 3 /s Mdm ) 


.(34) 


This is much less than the energy density of the CMB, 

Ecmb - 0.25 X 1012 f eV/cm3 , (35) 

implying that the distortions are below detection thresh¬ 
olds. We also note that although during matter dom¬ 
ination, the energy from DM annihilation is growing 
faster than the CMB energy density, the exponent of 9/2 
changes to 4 during radiation domination, and so the in¬ 
jected energy density from DM annihilation never is a 
substantial fraction of the CMB energy density. 

Could photons from DM annihilations, injected into 
an optically thin part of the spectrum, be detected 
as part of the diffuse photon background today? The 
transparency region from 10® to 10i° eV in Fig. ^ im¬ 
plies that these photons would have energies between 
10^ and lOi” eV today. Furthermore, Eq. for stan¬ 
dard parameters gives a present-day photon energy den¬ 
sity of E ^ 2 X 10-12 eV/cm®, implying a flux of 
0.5 X 10-2 eV/cm2/s/sr. Comparin g th is to the observed 
flux of 1 X 10® eV/cm2/s/sr we And that the 

fraction possibly due to DM annihilation is considerably 
below the uncertainties in the measurement. 

The residual ionization after recombination serves as a 
catalyst for the production of molecular hydrogen, first 
via 2H + ^ H 2 + at z ^ 500, and then by 

2H + e~ ^ H 2 + e~ at z 100. Molecular hydro¬ 
gen is important since it serves as a coolant, allowing for 
the collapse of the first cosmological objects. A higher 
ionization fraction could, in principle, create a greater 
abundance of molecular hydrogen, triggering collapse at 
an earlier epoch. The equation for the evolution of the 
molecular hydrogen fraction, / = n[H 2 ]/n[H\, is 

H{a)-^^^=keff{l-x-2f)nx, (36) 

where x is the ionization fraction and n = n[H]-\-n[E[^] + 
n[H 2 ] is the density of hydrogen atoms, and keff is 
the effective rate coefficient, approximately constant for 
X, / <C 1. For X <C 1, Eq. 1211 describes a catalyst-starved 
reaction; increasing x proportionally increases the H 2 
fraction. The standard recombination scenario produces 
an E [2 fraction of ~ IQ-® ^3 However, the fraction re¬ 
quired to efficiently cool halos is approximately IQ-®; pri¬ 
mordial H 2 , even with an enhanced ionization fraction, is 
too small to significantly alter structure formation. The 
required E 12 is produced in regions of high density, where 
the reaction rates are significantly higher. Unfortunately, 
this process is relatively insensitive to initial conditions 
|39| , making it difficult to constrain the H 2 and ionization 
fractions. 


VI. DISCUSSION 

We have considered the effect that DM annihilation has 
on the recombination history of the universe, and there¬ 
fore, on the CMB temperature and polarization power 
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spectra. We argued that the effect of DM annihilation 
on the IGM is well approximated by an injection of a 
fraction / of the rest mass energy of the DM particles 
into the IGM, where it is instantaneously used to heat 
and ionize the IGM. This “on the spot” approximation 
allowed us to compute the altered recombination history; 
the epoch of recombination is unchanged, but the residual 
ionization fraction increases. This broadens the visibility 
function, suppressing the temperature power spectrum, 
but enhancing the polarization power spectrum. Fur¬ 
thermore, the thicker visibility function shifts the peaks 
in TE and EE power spectra relative to the peaks in the 
TT power spectrum. 

Given the modified power spectra, we can ask whether 
the changes are detectable, or if they are degenerate with 
other cosmological parameters? A cosmic variance lim¬ 
ited survey is sensitive to an energy injection rate of 
Edm.o ~ 10“^"^ eV/s, probing masses ^100 GeV. Fur¬ 
thermore, these limits are attainable by CMB experi¬ 
ments designed to detect the polarization created by the 
stochastic gravitational wave background. The limits for 
Planck are about an order of magnitude worse, and the 
limits are significantly degraded for WMAP , due to its 
polarization sensitivity. However, WMAP will be able to 
constrain low mass (0(10 MeV — 1 GeV)) DM particles 
such as those proposed by 

We have kept our analysis as generic and idealized as 
possible, to ensure that our results are independent of 
the particulars of any DM model. We now consider some 
of the issues ignored by the analysis above. 

• How does one calculate / ? In order for a DM model 
to make a falsifiable prediction, it is important to 
be able to calculate / for a given model. A simple 
algorithm to compute / given all the decay channels 
and their branching ratios (admittedly a tedious 
task!) is : 

— Compute the and photon energy spectra 
resulting from an annihilation. 

— Evolve the photon and spectra with the 
processes in Eig. |21 

— Redshift the spectra to the next time and re¬ 
peat. 

One can obtain a qualitative picture of the re¬ 
sults by considering Fig. ^ For photons below 
^ 10® eV and electrons below 10^® eV, the en¬ 
ergy is efficiently deposited (~ 0.1 — 1) into 

the IGM. At higher energies, the efficiency is dic¬ 
tated by efficiency of pair production; however, 
even in this case, one would expect a non-negligible 
(~ 0.01 — 0.1) fraction to be converted into ioniza¬ 
tions and heat. One consequence is that one would 
expect / to generically increase with decreasing 
mass, weakening our ability to constrain the high¬ 
est masses but strengthening constraints on lower 
masses. Finally, to compare to plots such as Fig.[7| 
one must also include variations in {aAv) in /. 


While {gav) is only logarithmically dependent on 
Mdm for a thermal relic density of a single species, 
it can vary significantly if one includes the possi¬ 
bilities of multiple DM species or co-annihilations. 
Just how generic these processes are will depend 
on the class of theories being considered; initial at¬ 
tempts to answer such questions within the frame¬ 
work of the MSSM are in . 

• Is the standard recombination calculation accurate 
enough? Detecting DM annihilation by its effect on 
the recombination history requires that we under¬ 
stand the fiducial recombination physics to better 
than 1%. The current recombination calculation 
is accurate to ^ 1% ^3) although there recently 
have been claims of corrections due to two pho¬ 
ton processes ^3- Although the current calculation 
may not be accurate enough, we do not know of any 
theoretical limit that would prevent reaching the 
required accuracy. We should also emphasize that 
this accuracy is only necessary for detecting injec¬ 
tion energies ~ 10“^® eV/s; the current calculation 
is accurate enough to detect higher energies. 

• What about real-world complications? The limits 
on DM annihilation projected in Sec. ESI assumed 
an idealized CMB experiment with full sky cover¬ 
age, and no contamination to the primary CMB 
due to Galactic and extragalactic foregrounds. We 
briefly discuss these below; we however emphasize 
that these are not complications peculiar to mea¬ 
suring DM annihilation, but will affect any high 
precision polarization experiment (eg. detecting 
the gravitational wave background in the CMB). 

Decreasing the sky coverage reduces the number of 
available modes, increasing the errors due to cos¬ 
mic variance and degrading our ability to distin¬ 
guish between models. A full sky CMB experi¬ 
ment at these sensitivities might survey an effec¬ 
tive sky fraction of 50% due to Galactic and point 
source cuts, increasing the errors by -s/S; the impact 
of which is shown in Fig. El Even after exclud¬ 
ing the most contaminated regions of the mask, it 
will be necessary to separate Galactic and extra- 
galactic foregrounds using their frequency depen¬ 
dence. This problem has been examined in de¬ 
tail by a number of authors [eg. who demon¬ 
strate that this separation is possible on large angu¬ 
lar scales. Furthermore, extragalactic foregrounds 
start to dominate only on small angular scales 
[I ~ 2500), while the signal from DM annihilation 
is maximal on larger scales I < 1500 (Fig.EJ, sug¬ 
gesting that foregrounds are a tractable problem. 

We conclude by reiterating the complementarity be¬ 
tween the different probes of DM. As seen above, the 
GMB is able to put strong constraints on the cosmo¬ 
logical abundance of a light (<C 100 GeV) DM parti¬ 
cle that evades accelerator searches. At higher masses 
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(> 100 GeV), the injected energy from annihilations de¬ 
creases and the CMB is no longer competitive with other 
searches. However, given a DM candidate with a known 
mass and annihilation cross section, the injection energy 
constraints from the CMB translate into a constraint on 
the density of these particles. Importantly, this does not 
assume that the particle is a thermal relic, allowing us 
to constrain the particle content of the DM, including 
particles that produced via non-equilibrium processes. 

Understanding the properties of the dark matter re¬ 
mains one of the most important problems in cosmology 
today. And astrophysical probes, such as the one dis¬ 
cussed here, have an important role to play, both in de¬ 
tecting possible candidates and understanding their cos¬ 
mological abundance. 
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